fS^-nos! 


AFAL-TR-7g-16 
Volume  II 


FIRE  CONTROL  SYSTEM  ANALYSIS 

Volume  II  - Computer  Programming  Tasks 


Prepared  By: 

University  of  Dayton 
Research  Institute 
Dayton,  Ohio  45469 


November  1977 


Technical  Report  AFAL-TR-7^-16 , Vol  II 
Final  Report  for  Period  1 Nov  76  - 30  Sep  77 


Approved  for  public  release;  distribution  unlimited. 


AIR  FORCE  AVIONICS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

AI?  FORCE  SYSTEMS  COMMAND 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


D ^ 


NOTICE 


When  Gdveriment  drawings,  specifications , or  other  data  are  used  for  any  pur- 
pose other  than  in  connection  with  a definitely  related  Government  procurement 
operation,  the  United  States  Government  thereby  incurs  no  responsibility  nor  any 
obligation  whatsoever;  and  the  fact  that  the  government  may  have  formulated , 
furnished,  or  in  any  way  supplied  the  said  drawings,  specifications , or  other 
data,  is  not  to  be  regarded  by  implication  or  otherwise  as  in  any  manner  licen- 
sing the  holder  or  any  other  person  or  corporation,  or  conveying  any  rights  or 
permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any 
way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS) . At  NTIS,  it  will  be  avail- 
able to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


FOR  THE  COMMANDER 


( 


7 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list, 
or  if  the  addressee  is  no  longer  employed  by  your  organization  please  notify 
'/I A-’ ' 7~ / W-RAFB , OH  45433  to  help  us  maintain  a current  mailing  list". 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  se- 
curity considerations , contractual  obligations , or  notice  on  a specific  document. 


UNCLASSIFIED 


security  classification  of  this  page  flOian  Dmtt  Enttttd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  report  number  p.  GOVT  ACCESSION  NO. 

AFAL-TR-75-16,  Volume  II  j 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  title  (and  Su6l((l»i 

Fire  Control  System  Analysis 

Volume  II  - Computer  Programming  Tasks 

s.  TYPE  OF  report  4 PERIOD  COVERED 

Final  Report 

1 1 /1/76-9/30/77 

6.  PERFORMING  ORC.  REPORT  NUMBER 

7.  AUTHORfa) 

Carl  King 

6-  CONTRACT  OR  GRANT  NUHBCRfa) 

F33615-77-C-1056 

s 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

University  of  Dayton  Research  Institute 
Dayton,  Ohio  45469 

to.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  a WORK  UNIT  NUMBERS 

636O5F/T629O359 

n.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Avionics  Laboratory 
Wright-Patterson  AFB,  Ohio  45433 

12.  REPORT  DATE 

November  1977 

13.  NUMBER  OF  PAGES 

164 

14.  MONITORING  AGENCY  NAME  & ADDRESSfi/  dltlercnl  trom  Conlrotllng  Oflice) 

is.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

15a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  Cof  RepofO 

This,  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 

[ *7.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  from  Report) 

IB.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continue  on  reveree  aide  if  neceaaary  and  idenf//y  by  block  number) 

Kalman  Filter,  Target  State  Measurement,  Radar  Lag,  Multiple 
Digital  Scan  Converter,  Air  Combat  Evaluation  Algorithm 

20.  ABSTRACT  (Continue  on  reverae  aide  if  neceasary  and  Identify  by  block  number) 

The  following  three  tasks  are  included  in  Volume  I.  Gaussian  noise 
added  to  true  range  generated  by  various  kinds  of  simulated  tra- 
jectories served  as  "measurements"  which  were  input  to  several  kind 
of  Kalman  filters.  Filter  output  is  compared  graphically  and  by 
rms  deviations.  Independent  methods  were  devised  for  evaluating 
the  accuracy  of  target  state  vectors  (position,  velocity,  and 
acceleration)  obtained  from  a director  fire  control  system. 

(Conjinuedl 

DD  1473  EDITION  OF  t NOV  65  IS  OBSOLETE  UNCLASSIFIED 


security  classification  OF  THIS  PAGE  fH'han  Oaf*  Entarad) 


1 


J 


security  classification  of  This  PAGC(l*7i»n  0«l«  Enltrmd) 


UNCLASSIF. 


SECURITY  CLASSIFICATION  OF  THIS  PACEfHTi.n  Dmim  Enimfd) 


1 


TABLE  OF  CONTENTS 


SECTION 

1 

2 


3 


BACKGROUND 
TASK  DEFINITION 

2.1  REHOST  MODULAR  DITIGAL  SCAN  CONVERTER 
(MDSC)  CROSS  ASSEMBLER 

2.2  DEVELOP  AN  AERIAL  COMBAT  EVALUATOR 
(ACE)  ALGORITHM 

2.3  FIRE  CONTROL  DOCUMENTATION 
PROGRAM  SUPPORT 

3.1  REHOST  MODULAR  DIGITAL  SCAN  CONVERTER 
(MDSC)  CROSS  ASSEMBLER 

3.1.1  APPROACH 

3.1.2  VERIFICATION 

3.1.3  REHOST  ON  THE  CDC 

3.1.4  REHOST  ON  THE  PDP-11/55 

3.1.5  REFERENCES 

3.2  AIR  COMBAT  EVALUATOR  ALGORITHM 


3.2.1 

TASK  DEFINITION 

3.2.2 

ALGORITHM  DEVELOPMENT 

3.2.3 

THE 

RELATIVE  MOTION 

VECTOR 

3.2.4 

THE 

ANGLE  BETWEEN  THE  L AND 

R 

mo 

3.2.5 

THE 

MISS  DISTANCE 

3.2.6 

THE 

EXPECTED  NUMBER 

OF  HITS 

3.2.7 

THE 

ERROR  FUNCTION 

PAGE 

1 

3 

3 

3 

3 

5 

5 

5 

7 

7 

7 

10 

11 

11 

11 

15 

15 

16 
16 
16 


3.2.8  RESULTS 


19 


TABLE  OF  CONTENTS  (CONTINUED) 


SECTION  PAGE 

3 3.2.9  FLOWCHARTS  19 

3.2.10  INPUT-OUTPUT  19 

3.2.11  ALGORITHM  LISTING  28 

3.2.12  REFERENCES  28 

3.3  FIRE  CONTROL  DOCUMENTATION  29 

3.3.1  TASK  DEFINITION  29 

3.3.2  BASIC  SUPPORT  EFFORTS 

COMPLETED  30 

3.3.3  DOCUMENTATION  SUPPORT  30 

APPENDIX  A VERIFICATION  PROGRAMS  31 

APPENDIX  B COMPILED  PDP-11/45  LISTING  OF  CROSS 

ASSEMBLER  91 

APPENDIX  C ACE  ALGORITHM  145 

APPENDIX  D LAMARS  SUPPORT  PROGRAMS  149 


1 


LIST  OF  ILLUSTRATIONS 


FIGURE  PAGE 


1 

Approach  to  Rehosting 

6 

2 

Geometry  for  ACE  Algorithm 

12 

3 

Bullet  Stream  and  Target  Distribution 
Functions 

14 

4 

Miss  Distance  Diagram 

17 

5 

ERF(x)  and  the  Approximation 

to  ERF(x) 

20 

6 

Difference  Between  ERF  and  the  Approximation 
of  ERF 

21 

7 

Percent  Error  Between  ERF(x) 
Approximation  of  ERF(x) 

and  the 

22 

8 

Comparison  of  Miss  Distance  and  Expected 

Hits  for  Several  Relative  Motion  Vectors 

23 

9 

Descriptive  Flow  Chart 

24 

10 

ERF  Flow  Chart 

26 

SECTION  1 
BACKGROUND 


The  recognition  by  the  Air  Force  Avionics  Laboratory  (AFAL) 
of  the  need  for  improved  sensors  and  techniques  for  implementing 
fire  control  director  algorithms  has  resulted  largely  from  the 
conclusions  of  two  programs  - EXPO  V and  Fire/Fly.  EXPO  V 
is  the  latest  in  a series  of  man-in-the-loop  simulation  studies 
of  new  fire  control  and  gunsight  concepts.  The  Fire/Fly 
program  is  the  integrated  Fire  Control/Flight  Control  study 
which  is  investigating  methods  for  integrating  the  fire  control 
system  with  the  flight  control  system  to  improve  effectiveness 
while  increasing  survivability.  Both  programs  have  cited  the 
necessity  for  improved  sensors  as  well  as  improved  director 
algorithms . 

The  purpose  of  the  Fire  Control  Systems  analysis  effort, 
summarized  in  the  following  report,  was  to  examine  existing  fire 
control  systems  test  data  (Sight  Eval) , identify  error  sources, 
evaluate  and  modify  fire  control  algorithms,  and  perform 
programming  for  simulation,  weapon  system  investigation,  and 
validation . 

The  tasks  contained  in  this  volume  were  for  computer 
programming  support. 
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SECTION  2 
TASK  DEFINITION 

The  broad  categories  of  effort  contained  in  the  Statement  ] 

1 

of  Work  were  refined  by  coordination  with  Captain  J.  Silverthorn,  i 

RWT-2,  and  resulted  in  the  following  series  of  tasks.  j 

2.1  REHOST  MODULAR  DIGITAL  SCAN  CONVERTER  (MDSC)  CROSS  H 

ASSEMBLER  ' j 

i 

The  Hughes  Aircraft  Company  Modular  Digital  Scan  Converter  ] 

(MDSC)  forms  an  important  part  of  the  director  gunsight  flight  test. 

It  will  be  performing  important  computations  and  Heads  Up  Display 
(HUD)  symbol  generation.  It  is  imperative  that  a capability  exist 
at  Tyndall  Air  Force  Base  (TAFB)  to  modify  the  MDSC  assembly 
language  program  in  order  to  change  computations  or  symbols. 

To  accomplish  this,  an  existing  MDSC  cross  assembler  will  be 

rehosted  on  either  a PDP-11  or  ROLM  16/64  computer,  or  both.  ! 

i 

2.2  DEVELOP  AN  AERIAL  COMBAT  EVALUATOR  (ACE)  ALGORITHM  j 

An  Aerial  Combat  Evaluator  (ACE)  algorithm  was  used  on  j 

the  Sight  Eval  program  and  shown  to  be  valuable.  It  indicates  ‘ 

! 

to  the  pilot  the  expected  number  of  bullets  that  would  have  hit  j 

the  target  had  he  fired  the  gun.  As  a result,  it  provides  | 

feedback  even  during  "dry  run"  passes.  The  algorithm  used  in 
Sight  Eval  and  the  algorithm  developed  in  a previous  contract  ; 

will  be  compared,  the  best  selected,  and  then  implemented  on 

the  ROLM  computer . 1 

2.3  FIRE  CONTROL  DOCUMENTATION 

A significant  amount  of  documentation  of  previously 
developed  algorithms  needs  to  be  performed.  Many  of  these 
algorithms  are  being  used  in  the  director  flight  test. 
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SECTION  3 
PROGRAM  SUPPORT 

3.1  REHOST  MODULAR  DIGITAL  SCAN  CONVERTER  (MDSC) 

CROSS  ASSEMBLER 

The  Modular  Digital  Scan  Converter  (MDSC)  computer  was 
developed  by  Hughes  Aircraft  Company  for  use  in  airborne  fire 
control  systems.  Since  the  MDSC  computer  does  not  have  an 
assembly  language  compiler,  programming  is  done  through  a 
cross  assembler.  (A  cross  assembler  is  a program  that 
generates  machine  code  for  a given  computer  from  the  assembly 
language  of  another  computer.)  The  cross  assembler  for  the  MDSC 
computer  was  initially  hosted  on  the  Sigma  V computer  at  f 

Hughes  and  all  programming  had  to  be  done  through  this  system. 

To  expedite  work  on  programming  on  the  MDSC  computer,  the  ARAL 
of  Wright-Patterson  Air  Force  Base  (WPAFB)  contracted  with  the 
University  of  Dayton  Research  Institute  (UDRI)  to  rehost  the 

MDSC  Cross  Assembler  to  the  PDP-11  computer  and  the  ROLM  16/64  i 

computer.  These  computers  are  on-site  at  TAFB  where  the  Cross  | 

Assembler  program  is  to  be  implemented. 

3.1.1  Approach 

The  rehosting  process  is  illustrated  in  Figure  1.  The 
first  step  in  rehosting  the  FORTRAN  coded  MDSC  Cross  Assembler  to 

the  PDP-11  and  ROLM  16/64  computers  was  to  rehost  it  to  the  I 

CDC  computer  at  WPAFB.  The  step  of  rehosting  the  MDSC  Cross 
Assembler  to  the  CDC  computer  system  was  taken  because  of  the 
familiarity  and  availability  of  this  system  to  the  UDRI  personnel. 

An  added  advantage  of  this  step  was  that  the  CDC  system  checked 

the  operation  of  all  subsections  of  a program.  Once  the  MDSC  1 

Cross  Assembler  was  operational  on  the  CDC  system,  work  was 

directed  to  the  task  of  rehosting  it  on  the  PDP-11  and  ROLM 

16/64  computers. 

The  PDP-11/55  computer  available  at  TAFB  has  the 
RSX-llM  operating  system  and  a FORTRAN  IV-PLUS  compiler.  UDRI 
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Figure  1.  Approach  to  Rehosting 


was  unable  to  locate  a PDP-11/55  computer  at  WPAFB,  but  a 
PDP-11/45  with  the  same  operating  system  and  compiler  was 
located.  This  proved  to  be  a satisfactory  substitute  since 
the  two  systems  are  operationally  the  same.  The  rehosting  of 
the  Cross  Assembler  to  the  PDP-11/45  was  completed  in  August 
1977  and  sent  to  TAFB  for  implementation. 

Although  the  initial  plans  were  to  rehost  the  Cross 
Assembler  to  both  the  PDP-11  and  the  ROLM  16/64  computers, 
only  the  rehosting  to  the  PDP-11  was  completed.  Rehosting 
to  the  ROLM  16/64  computer  could  not  be  initiated  because  of 
scheduling  problems  on  this  computer.  The  Air  Force  Flight 
Dynamics  Laboratory  (AFFDL)  had  priority  on  this  computer  and 
its  current  programs  consumed  practically  all  of  the  avail- 
able time. 

3.1.2  Verification 

Verification  of  the  rehosted  MDSC  Cross  Assembler  was 
accomplished  by  running  five  documented  assembly  language 
programs  on  each  of  the  versions  and  comparing  the  output  values 
to  the  verification  values.  Two  of  the  verification  programs 
were  written  by  Hughes  (Reference  1)  and  three  by  AFAL 
(Reference  2) . Four  of  the  five  agreed  perfectly  with  the 
verification  listings.  The  output  of  one  of  the  AFAL  programs 
did  not  agree  fully  with  the  expected  output.  However,  this 
was  not  considered  a serious  problem  since  this  program  has 
not  been  completely  debugged.  The  listing  of  these  five 
verification  programs  as  compiled  by  the  PDP-11  version  of  the 
MDSC  Cross  Assembler  are  shown  in  Appendix  A. 

3.1.3  Rehost  on  the  CPC 

Rehosting  the  MDSC  Cross  Assembler  program  from  the 
Sigma  V computer  to  the  CDC  computer  system  required  several 
modifications  because  of  hardware  differences.  The  major 
difference  in  the  two  computer  systems  is  word  size.  The 
Sigma  V uses  a 32-bit  word  with  four  bytes  per  word  and  the 
CDC  uses  a 60-bit  word  with  10  bytes  per  word.  To  circumvent 
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this  problem,  Sigma  V coding  was  changed  to  the  CDC  format  of 
10  bytes/word. 

Another  difference  in  the  two  systems  is  that  the 
Sigma  V is  a two's  complement  computer  and  the  CDC  is  a one's 
complement  computer.^  It  was  found  that  this  problem  could  be 
overcome  by  modifying  the  masking  subprogram  to  make  the  CDC 
function  as  a two's  complement  computer. 

Modifications  that  had  to  be  made  to  further  make  the 
Cross  Assembler  compatible  to  the  CDC  system  were  as  follows. 

(1)  A program  specification  statment  had 
to  be  added  as  the  first  statement  of  the 
main  program. 

(2)  Variable  and  subroutine  names  had  to 
be  reduced  to  seven  or  less  characters. 

(3)  The  end-of-file  check  and  READ  statement 
had  to  be  modified. 

(4)  Logical  unit  numbers  for  input  and 
output  had  to  be  less  than  100. 

(5)  All  DO  statements  had  to  be  changed 
so  that: 

a)  no  computations  occurred  within  the 
DO  statement. 

b)  all  index  increments  were  positive, 

c)  the  index  started  at  an  integer 
greater  than  or  equal  to  1. 

d)  the  branch  statements  back  into  a 
DO  loop  had  to  be  eliminated. 

(6)  External  function  names  had  to  be 
changed  to  agree  with  those  used  in  the 
CDC  FORTRAN  compiler. 

(7)  FORMAT  statements  were  changed  to  make 
the  output  more  readable. 


A one's  complement  computer  represents  negative  integers  as  the 
complement  bit-by-bit  of  the  positive  integers.  A two's 
complement  computer  represents  negative  integers  as  a one's 
complement  with  the  number  one  added  to  the  negative  integer. 
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Several  definite  coding  errors  were  detected  in  the 
Cross  Assembler  program  and  corrected  in  the  rehosting  process. 
Several  of  the  major  errors'” occurred  in  the  assembly  language 
symbolic  names  and  the  machine  language  for  each.  These 
errors  were  corrected  to  agree  with  Reference  1.  In  addition, 
the  tabbing  subprogram  code  was  changed  to  tab  ove . commas 
and  blanks  in  the  assembly  language  cards.  As  a result  of 
rehosting  the  Cross  Assembler  program  to  the  CDC,  the  operation 
of  each  subprogram  was  checked  and  demonstrated  to  be  accurate. 

3.1.4  Rehost  on  the  PDP-11/55 

The  major  changes  required  to  rehost  the  MDSC  Cross 
Assembler  to  the  PDl -11/55  computer  were  similar  to  those 
required  for  the  rehost  to  the  CDC  computer.  Machine  com- 
patibility was  not  a problem  in  i.his  rehosting  since  the  PDP-11/55 
computer  with  the  RSX-llM  operating  system  and  FORTRAN  IV-PLUS 
compiler  mimics  the  word  structure  of  the  Sigma  V computer. 

In  rehosting  the  CDC  compatible  Cross  Assembler  to 
the  PDP-11,  the  following  changes  were  made. 

(1)  Variable  and  subroutine  nciTues  had  to  be 
reduced  to  six  or  less  characters. 

(2)  DATA  statements  had  to  be  altered  to 
eliminate  implied  DO  loops. 

(3)  The  end-of-file  check  in  the  READ 
statement  had  to  be  modified. 

(4)  The  FORMAT  had  to  be  changed  to  output 
a PDP-11  compatible  output. 

(5)  External  function  names  had  to  be 
changed  to  the  corresponding  names  in  the 
PDP  compiler. 

(6)  Multiple  assignment  statments  had  to 

be  recoded  so  that  there  was  no  more  than  one 
equal  sign  in  each  line  of  code. 

(7)  A routine  had  to  be  coded  to  handle 
the  decoding  of  hexidecimal  input  integers 
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as  characters  to  the  respective  numerical 
integer  values  used  to  compute  machine 
language  instruction  from  assembly  language 
input. 

(8)  DOUBLE  PRECISION  statements  removed 
in  rehosting  to  the  CDC  system  had  to  be 
replaced . 

In  compiling  the  Cross  Assembler  on  the  PDP-11 
computer,  it  was  found  that  the  Cross  Assembler  program  exceeded 
the  capacity  of  the  machine.  The  maximum  capacity  of  the  PDP-11 
is  32  K words  of  16-bit  length.  To  obviate  this  problem,  the 
concordance  table  of  assembly  language  labels  was  deleted  from 
the  program.  After  the  above  modifications  were  made  and  the 
program  shown  to  be  operational,  punched  card  decks  of  the 
Cross  Assembler  and  its  verification  programs  were  delivered 
to  AFAL.  Appendix  B is  a listing  of  the  MDSC  Cross  Assembler 
program. 

3.1.5  References 


(1)  MDSC  Programmer's  Reference  Manual  Volume  4,  Hughps 

Aircraft  Company,  Culver  City,  California,  HAC 
Reference  no.  D2385,  March  1976. 

(2)  AN/UYK-30  Reference  Manual,  Hughes  Aircraft  Company, 

Culver  City,  California,  Report  no.  P76-148, 
Ref.  A,  November  1976. 


3.2 


AIR  COMBAT  EVALUATION  ALGORITHM 


3.2.1  Task  Definition 

The  University  of  Dayton  Research  Institute  has  developed 
an  Air  Combat  Evaluation  (ACE)  algorithm  to  simulate  the  number 
of  hits  scored  during  air  encounters  between  an  attacking  air- 
craft and  a target  aircraft.  This  algorithm,  when  integrated 
into  the  flight  control  system  of  an  attack  aircraft,  will  allow 
the  pilot  to  experience  combat  conditions  involving  evasive  maneu- 
vers without  using  live  ammunition.  The  obvious  consequence  of 
the  use  of  this  algorithm  is  improved  aircraft-aircraft  combat 
training  techniques. 

The  ACE  algorithm  uses  a statistical  approach  to  compute 
the  expected  number  of  hits  scored  by  the  attacking  aircraft  on 
the  target  aircraft.  This  information  is  recorded  and  available 
to  the  pilot  on  a real-time  basis.  The  position  of  the  target 
aircraft  is  obtained  from  the  fire  control  system  of  the  attacking 
aircraft.  Specific  inputs  to  the  ACE  algorithm  are  azimuth  and 
elevation  of  the  computed  projectile  stream  at  target  range  with 
respect  to  the  target. 

3.2.2  Algorithm  Development 

The  ACE  algorithm  is  a modification  of  the  Uniform  Normal 
Algorithm  (Appendix  1;  Ref.  1).  ACE  is  coded  in  FORTRAN  as  a 
subroutine  for  implementation  on  the  ROLM  16/64  computer.  However, 
because  of  the  unavailability  of  the  ROLM  system  the  major  coding 
effort  was  done  using  the  CDC  computer  system. 

The  geometry  for  the  ACE  algorithm  is  shown  in  Figure  2 . 

The  target  plane  is  normal  to  the  line  of  sight  from  the  attack 
aircraft  to  the  target  with  its  origin  at  the  center  of  the  tar- 
get. The  point  of  intersection  of  the  previous  projectile  stream 
(one  computer  cycle  before)  with  the  target  plane  is  represented 

by  the  vector  R and  the  point  of  intersection  of  the  current 
mo 
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Figure  2 . Geometry  for  ACE  Algorithm 
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projectile  stream  by  R^.  The  relative  motion  of  the  projectile 
stream  is  represented  by  L.  The  left-handed  coordinate  system 
is  chosen  with  the  y-axis  parallel  to  L and  positive  x down. 

The  trajectory  of  the  projectile  stream  is  that  of  a 
nominal  projectile.  The  dispersion  o'  projectiles  about  the 
projectile  stream  is  assumed  to  be  uniform  along  L and  random 
and  normally  distributed  in  the  target  plane  as  shown  in  Figure 
2.  The  density  function  for  the  projectile  stream  is 


e~  20  2 
Lo  2tt 

O 


= O 


|yl<  L/2 


|y| > L/2 


(1) 


where  is  the  standard  deviation  of  this  projectile  stream,  in 

£3 

radians,  and  L is  the  magnitude,  in  radians,  of  the  relative 
motion  vector  L. 

The  target  is  represented  by  a bivariate  normal  distri- 
bution in  the  target  plane  as  shown  in  Figure  3 . The  density 
function  for  the  target  is 


(x-x^) 2+(y-y^) 2 


f(x,y)  = 


(2] 


where  o .j,  is  the  standard  deviation  of  the  target,  in  feet,  and 


T' 


y,  in  feet.  The  standard  deviation  of  the  target  in 


radians  is  then  given  by 


TR 


Using  this  notation  the  density  function  of  the  target  can  be 

2 2 

written  (x-x^)  +(y-y^) 

2o__2 

® (3) 

where  x,  x^,  y,  y^  are  now  also  radians. 
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The  expected  number  of  hits  on  the  target  is  given  by  the 
product  of  the  number  of  projectiles  in  the  interval  L and  the 
probability  that  one  projectile  will  hit  the  target.  The 
probability  that  one  projectile  will  hit  the  target  is  the 
intersection  of  the  two  density  functions.  Therefore,  the 
expected  number  of  hits  is 

^ L/o 


E = N 


B 


If  f,p(x,y)fgg(x,y)dy  dx 


(4) 


-00  -L/, 


Preliminary  to  computing  expected  number  of  hits,  several 
important  intermediate  values  must  be  found.  The  principle 
computations  are  described  below. 


3.2.3 


The  Relative  Motion  Vector 


The  relative  motion  vector  defines  the  motion  of 
the  projectile  stream  between  the  previous  position  and  the 
current  position.  The  magnitude,  L,  of  this  vector  is  the  dis- 
tance between  successive  projectile  streams  at  target  range 
measured  relative  to  the  target  position  in  the  plane.  Thus, 

L is  given  by 


L = R - R 
mo  m 


(5) 


where  R is  the  previous  position  vector  and  R is  the  current 
mo  i.  IT 


position  vector. 


3.2.4 


The  Angle  Between  the  L and  R 


mo 


The  angle  between  the  vectors  L and  R is  repre- 
^ mo 


sented  by  6 in  Figure  1 and  is  computed  by  the  dot  product  rela- 
tionship 

L‘R 


cos  6 = -rt^ 


mo 


(6) 


mo 


T1 
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1 

j 

3.2.5  The  Miss  Distance 

The  miss  distance,  with  components  and  y^,  is  | 

the  line  in  the  target  plane  from  the  target  to  the  midpoint  of  j 

L,  as  shown  in  Figure  4.  The  component  is  the  normal  distance  1 

from  the  target  to  L and  is  given  by  j 

] 

^T  " j 


The  component  y^  is  along  L and  is  given  by 


3.2.6  The  Expected  Number  of  Hits 


The  expected  number  of  hits  is  calculated  from 
Equation  (4)  which  has  been  integrated  to  obtain  the  following 
closed  form  expression 


where  Ng  is  the  number  of  projectiles  at  target  range  per  com- 
puter cycle  time  (fire  rate  x cycle  time) , a is  the  angular 
standard  deviation  of  the  target  (radians) , a„  is  the  angular 
standard  deviation  of  the  projectile  stream  (radians) , L is  the 
magnitude  of  the  relative  motion  vector  (radians) . 

3.2.7  The  Error  Function 


The  Error  Function  used  in  the  ACE  algorithm  is  given  by 


o 


du 


(10) 


\ 
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Figure  4.  Miss  Distance  Diagram 
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This  function  has  the  following  properties: 


ERF  (-X)  = -ERF(x) 

ERF  (o°)  = 1/2 
ERF(O)  = 0 
ERF  (x>3.5)«l/2 

The  error  function  described  in  (10)  cannot  be  evaluated  in  closed 
forM  and,  therefore,  its  value  must  be  found  numerically. 


Four  methods  of  approximating  values  for  the  error 
function  were  considered. 

(1)  The  use  of  a table  of  values  for  the  error  function 
at  specified  argument  values. 

(2)  .A  series  expansion  of  the  error  function. 

(3)  Numerical  integration  of  the  Equation  10. 

(4)  A least  square  polynomial  representation  of  the 
error  function. 


These  four  methods  of  approximating  values  of  the  error  function 
were  compared  relative  to  computer  space  and  time  requirements. 

As  a result  of  this  comparison  it  was  concluded  that  the  best 
method  for  this  subroutine  was  a least  square  polynomial  fit. 

A regression  technique  was  used  to  fit  third  through  sixth  degree 
polynomials  to  (10)  the  error  function.  A third  degree  poly- 
nomial was  chosen  for  the  fit  with  appropriate  linear  corrections 
for  small  and  large  x.  The  specific  error  function  approximation 
is 


2 3 

ERF(x)  = Kg+Kj^x+K2X 


(11) 


where. 


for  X e [0 . 0 , 0.2] 


= 0.39629855 


X £(0.2,  3.5) 


Kq  = K2  = K3  = 0 

Kq  = -0.01322336 
= 0.49613046 
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= -0.15942666 
= 0.01698544 

X e (3.5,  0°)  Kq  = 0.5 

Xl  = Kj  = K3  =0 

The  least  square  fit  of  linear  polynomials  for  0_<x<0.2  and 

x^3.5  was  done  to  improve  the  fit  in  these  regions.  The  error 

function  and  the  corresponding  approximated  values  are  plotted 

in  Figure  5.  The  actual  differences  between  the  true  values 

of  the  error  function  and  the  approximated  values  is  more  clearly 

shown  in  Figure  6.  The  corresponding  percent  difference  is 

shown  in  Figure  7 . The  improvement  in  the  approximation  in  the 

regions  0_<x<0.2  and  x^3 . 5 can  be  seen  in  these  figures.  ^ 

3.2.8  Results 

Results  obtained  from  the  ACE  algorithm  for  selected 
input  parameters  are  shown  in  Figure  8.  The  number  of  expected 
hits  is  plotted  versus  miss  distance  in  this  figure.  No 
empirical  data  is  available  at  this  time  for  comparison.  If 
there  is  a discrepancy  in  the  results  obtained  and  the  antici- 
pated results,  a detailed  study  of  the  approximation  to  the 
error  function  should  be  considered,  because  the  given  approxi- 
mation could  lead  to  rather  large  errors  in  the  expected 
number  of  hits. 

3.2.9  Flowcharts 

Figure  9 shows  the  descriptive  flow  chart  for  the  ACE 
algorithm  and  Figure  10  shows  the  flow  chart  for  the  error  func- 
tion . 

3.2.10  Input-Output 


INPUTS 

VARIABLE  NAME 

UNITS 

DESCRIPTION 

AZ 

Radians 

The  azimuth  position  of  the 
projectiles  at  target  range 
respect  to  the  target 

computed 

with 

EL 

Radians 

The  elevation  position  of  the  com- 

puted  projectiles  at  target  range 
with  respect  to  the  target 


I 

'I 


Figure  5.  ERF(x)  and  the  Approximation  to  ERF(x) 


Figure  7.  Percent  Error  Between  ERF(x)  and  the 
Approximation  to  ERF(x). 
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COMPUTE 
MISS  VECTOR 


DESCRIPTION 

The  range  of  the  target  with  respect 
to  the  attacker 

The  initialization  flag  that  either 
initializes  or  runs  the  subprogram. 
If  IFLAG  = 1,  then  the  program  is 
run  and  computes  the  expected  number 
of  hits. 

If  IFLAG  1,  then  the  program 
initializes  several  constants  to 
be  used  in  the  program. 

The  input  variables  are  in  common  and  the  calling  program 
must  have  a common  statement  as  follows  to  input  values  to  this 
subprogram; 

COJiMON/INPT/AZ,  EL,  RANGE,  IFLAG  . 

OUTPUT 


VARIABLE  NAME  UNITS 

RANGE  Feet 

IFLAG 


VARIABLE  NAME  DESCRIPTION 

EXPHTS  This  value  is  the  expected  number 

of  hits  of  the  projectiles  on  the 
target  when  both  are  at  target  range. 
It  is  the  expected  number  of  hits  per 
computer  cycle. 

This  output  variable  is  in  common  and  the  calling  program 
must  have  a common  statement  as  follows  to  output  values  from  this 
subprogram; 

COMMON/OUTPT/EXPHTS  . 


INPUT 

CONSTANTS 

VARIABLE  NAME 

UNITS 

DESCRIPTION 

BULN 

The  number  of  projectiles  at 
range  per  computer  cycle  time 
rate  x cycle  time) 

target 

(fire 

SIGT 

Feet 

The  standard  deviation  of  the 

target 

SIGBS 

Radians 

The  standard  deviation  of  the 
jectile  stream  or  projectile 
persion. 

pro- 

dis- 
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The  constant  variables  are  common,  and  the  calling 
program  must  have  a common  statement  as  follows  to  input  these 
constant  values  to  the  subprogram; 

COMMON/CONST/BULN,  SIGT,  SIGBS  . 

INTERNAL  CONSTANTS 


VARIABLE  NAME 


SQ2PI 


UNITS 


AZOLD 


ELOLD 


Radians 


Radians 


Radians' 


DESCRIPTION 
Square  root  of  2 x tt 
Exponential 

Previous  value  of  the  azimuth 
position  of  the  computed  projectiles 
at  target  range  with  respect  to 
the  target 

Previous  value  of  the  elevation 
position  of  the  computed  pro- 
jectile at  target  range  with 
respect  to  the  target 

SIGBS^ 


These  variables  are  common  only  to  this  subprogram.  The 
later  three  are  in  common  and  are  assigned  their  values  at  the 
initialization  stage  of  the  subprogram.  The  common  is  as  follows; 

C0MM0N/ACEC0M/SQ2PI,  E,  AZOLD,  ELOLD,  W2  . 

The  first  several  values  are  constant  at  all  times  and 
therefore  are  in  data  statements  as  follows; 

DATA  SQ2PI/2. 5066283/  , 

DATA  E/2.7182818/  . 

3.2.11  Algorithm  Listing 

The  listing  of  the  ACE  algorithm  as  a coded  FORTRAN  sub- 
routine is  contained  in  Appendix  C. 

3.2.12  References 


Edwards,  Verlan  E.  and  Ted  G.  Johnson,  Air-to-Air  Gunfire 
Control  System  Evaluator,  Honeywell,  Inc.,  Technical 
Report  AFAL-TR-73-20 , November  1972. 
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3.3 


FIRE  CONTROL  DOCUMENTATION 


3.3.1  Task  Definition 

The  Large  Amplitude  Aerospace  Research  Simulator  (LAMARS) 
facility,  which  is  supported  by  the  AFFDL  at  WPAFB,  is  designed 
to  simulate  the  total  environment  of  an  aircraft  in  flight. 

The  LAMARS  facility  consists  of  a cockpit  mounted  on  a large 
movable  arm,  various  television  picture  projectors  and  a 
computer  controlled  servo-system.  The  cockpit  area  can  be 
modified  to  meet  the  cockpit  specifications  of  the  aircraft 
being  simulated. 

Aircraft  motion  is  simulated  by  a combination  of  arm 
movement  and  motion  picture  projection  of  terrain.  The  move- 
ment of  the  arm  is  activated  through  the  controls  in  the  cock- 
pit as  the  pilot  "flys"  the  aircraft.  The  terrain  over  which 
the  pilot  is  flying  is  simulated  by  projecting  televised 
pictures  of  a terrain  onto  a screen  that  surrounds  the  cockpit 
area.  When  the  facility  is  being  used  to  simulate  fighter 
aircraft  in  combat,  the  image  of  a target  is  projected  onto  the 
screen.  The  flight  path  of  the  target  aircraft  is  input  to 
the  system. 

The  LAMARS  facility  was  modified  to  simulate  the  F-106 
fighter  aircraft.  The  major  aircraft  responses  and  terrain 
along  with  the  target  were  provided  by  AFFDL  and  the  fire 
control  algorithms  were  provided  by  AFAL.  The  integration  of 
the  AFAL  algorithms  took  place  on  the  ROLM  16/64  computer  on 
the  computer  deck  at  LAMARS  facility.  The  ROLM  16/64  computer 
was  used  because  this  computer  is  part  of  the  air-borne  computer 
system  in  the  F-106  aircraft.  The  ROLM  16/64  used  in  the  LAMARS 
facility,  is  connected  to  several  other  computers  on  the  com- 
puter deck  by  use  of  a Direct  Memory  Access  (DMA)  channel.  The 
inputs  to  the  fire  control  algorithms  are  passed  from  the  air- 
craft simulation  computer  through  the  DMA  to  the  ROLM  16/64. 

The  outputs  of  the  fire  control  algorithms  are  then  passed 
from  the  ROLM  16/64  through  the  DMA  to  the  symbol  generator 
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computer  to  input  the  symbology  to  the  Heads  Up  Display 
for  viewing  by  the  pilot. 


The  addition  of  the  fire  control  algorithms  to  the  simu- 
lation facility  allows  the  pilot  to  use  the  symbology  on  the 
HUD.  The  pilot  can,  as  in  an  actual  flight,  select  what  is  to 
be  displayed  on  the  HUD  and  then  use  the  symbology  to  track  the 
image  of  the  target  on  the  screen.  This  addition  to  the  LAMARS 
facility  adds  to  the  realism  of  flight  to  the  pilot  and 
helps  to  train  the  pilot  to  effectively  fly  a fighter  aircraft 
to  track  a target  while  in  flight. 

3.3.2  Basic  Support  Efforts  Completed 

The  role  of  the  University  of  Dayton  Research  Institute 
in  the  integration  of  the  fire  control  algorithms  was  basically 
one  of  a supportive  nature.  The  AFAL  planned  the  entire 
integration  package  using  the  most  common  top-down  programming 
approach.  The  UDRI  assisted  in  integrating  five  of  the  fire 
control  algorithms  into  the  planned  integration  set  by  the  AFAL. 
Each  algorithm,  coded  in  FORTRAN  for  the  ROLM  16/64  computer, 
had  to  be  adapted  for  the  integration.  The  adaptation  of  each 
algorithm  required  the  specification  of  the  inputs,  outputs, 
ana  constants  by  the  appropriate  common  blocks  as  set  by  AFAL. 

The  completed  algorithms  for  the  integration  are  listed 
in  Appendix  D.  These  listings  are  the  compiled  versions  of  the 
five  algorithms  adapted  for  the  LAMARS  integration. 

3.3.3  Documentation  Support 

Sample  documentation  was  provided  to  Captain  Silverthorn 
to  assist  in  documenting  Hot  Line  Gunsight  (HLGS)  and  Lead 
Computing  Optical  Sight  (LCOS) . 
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APPENDIX  A 

VERIFICATION  PROGRAMS 

Five  assembly  language  programs  were  obtained  from 
Reference  1 (see  Paragraph  3.1.5  of  this  volume)  to  verify 
the  operation  of  the  Cross  Assembler.  These  programs  were 
v/ritten  for  the  MDSC  computer  and  compiled  by  the  PDP  version 
of  the  Cross  Assembler  on  the  PDP-11/45  at  WPAFB.  Verification 
Programs  1 and  2 are  those  taken  from  Reference  1.  Verification 
Programs  3,  4,  and  5 are  those  supplied  by  the  AFAL  branch 
which  deals  with  the  MDSC  computer.  Verification  Program  3 
has  the  complete  test  of  the  assembly  language  instruction  set 
and  all  modes  of  each  instruction.  It  should  be  noted  that 
Verification  Program  4 has  an  error  indicated.  Also,  errors 
have  been  detected  in  the  verification  listing.  Consequently, 
there  is  a question  as  to  the  validity  of  this  program. 

Verification  Programs  1,  2,  and  5 are  typical  of  the  programs 
used  in  the  MDSC  computer.  Verification  Programs  1,  2,  3,  and 
5 have  been  checked  out  with  the  verification  listings. 

The  output  format  consists  of  4 columns. 

Column  1 is  the  card  or  line  number  of  the 
assembly  language  instruction 

Column  2 is  the  location  counter  (octal) 

Column  3 is  the  machine  code  (octal) 

Column  4 is  the  input  assembly  language 
instruction  which  generated  the  machine  code 

A listing  of  each  of  the  five  verification  programs  follows. 
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APPENDIX  B 

COMPILED  PDP-11/45  LISTING  OF  CROSS  ASSEMBLER 
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Following  is  a listing  of  the  MDSC  Cross  Assembler 
program  coded  in  FORTRAN  as  compiled  by  the  PDP-11/45  computer 
at  WPAFB  with  the  RSX-llM  operating  system  with  the  FORTRAN  IV- 
PLUS  compiler. 
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APPENDIX  C 
ACE  ALGORITHM 

Following  is  a listing  of  the  ACE  Algorithm  subprogram 
coded  in  FORTRAN  as  compiled  by  the  CDC  CYBER-74  series  computer 
system. 
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APPENDIX  D 

LAMARS  SUPPORT  PROGRAMS 

Following  are  the  listings  of  the  five  subpiograms 
as  compiled  by  the  CDC  CYBER-74  series  computer  system  for  the 
LAMARS  fire  control  integration  to  be  implemented  on  the  ROLM 
16/64  computer. 
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